Asymmetric exclusion model with impurities. 
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An integrable asymmetric exclusion process with impurities is formulated. The model displays 
the full spectrum of the stochastic asymmetric XXZ chain plus new levels. We derive the Bethe 
equations and calculate the spectral gap for the totally asymmetric diffusion at half filling. While 
the standard asymmetric exclusion process without impurities belongs to the KPZ universality class 
with a exponent § , our model has a scaling exponent § . 
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One-dimensional 3-state quantum Hamiltonians and 
master equations of reaction-diffusion processes have 
played an important role in describing strongly corre- 
lated electrons and non-equilibrium statistical mechanics 
in the last decades, mainly due to their intrinsic and non- 
trivial many-body behavior. Remarkably, in one dimen- 
sion several models in this category are exactly solvable, 
as the spin-1 Sutherland [TJ and t-J [5j models, and the 
asymmetric diffusion of two types of particles 3 . In its 
formulation in terms of particles with two global conser- 
vation laws, these models describe the dynamics of two 
types of particles on the lattice, where the total number 
of particles of each type is conserved separately. In order 
to ensure integrability, all known models in this class sat- 
isfies some particle-particle exchange symmetries [H [5] . 
Recently, we introduced a new class of 3-state model in 
the context of high energy physics that is integrable de- 
spite it do not have particle-particle exchange symmetry 
[5]- The quantum version is closely related to the strong 
regime of the t-U Hubbard model and the XXC model [B] . 
In this work we formulate a stochastic model related to [5] 
that describes an asymmetric exclusion process (ASEP) 
with impurities. Although our model can be solved by 
the coordinate Bethe ansatz, we are going to formulate a 
new matrix product ansatz (MPA) [H [7] due its simplic- 
ity and unifying implementation for arbitrary systems. 
This new MPA introduced in [4j [7] can be seen as a ma- 
trix product formulation of the coordinate Bethe ansatz 
and it is suited to describe all eigenstates of integrable 
models. We solve this model with periodic boundary 
condition through the MPA and we analyze the spectral 
gap for some special cases. Our model displays the full 
spectrum of the ASEP without impurities [8] plus new 
levels. The first excited state belongs to these new levels 
and has unusual scaling exponents. Although the ASEP 
without impurities belongs to the KPZ universality class 
[9] (dynamic exponent |), our new model displays a scal- 
ing exponent § . 

The model we propose describes the dynamics of two 
types of particles (type 1 and 2) on the lattice, where the 
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total number of particles of each type is conserved. In 
this model if the neighbors sites are empty, particles of 
type 1 can jump to the right or to the left with rate Tq ° 
and rj q , respectively. Particles of type 2 "impurities" 
do not jump to the neighbors sites if they are empty, but 
can change positions with neighbors particles of type 1 
with rates Tj \ and T 2 \ if the particle 1 is on the left or 
on the right, respectively. There are not particle-particle 
exchange symmetry because particles 1 and 2 have differ- 
ent kinds of dynamics. The stochastic model we consider 
on a ring of perimeter L is related to the Hamiltonian 
H = HW + #(2) where (a = 0, 2) 
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and E a 'P (a, ft = 0, 1, 2) are the usual 3x3 Weyl matrices 
with i,j elements {E l > m ), ,. = 5i,i5 m ,j, and Tif + Tf I = 
1 . Finally, the model (IT]) is related to the tracer diffusion 
model |3] by an interchange of labels and 1, however 
the exactly solutions of these models describe different 
sectors of particles occupation. 

The conservation of particles supplemented by the pe- 
riodic boundary condition of (flj) imply that the total 
number of particles n\ , n-i on class 1 and 2 as well the mo- 
mentum P = ^fl (J = 0, 1, . . . , L — 1) are good quantum 
numbers. The eigenstates I 1 ?) of the eigenvalue equation 
-ff|vP) = £ ni - n2 y\>) belonging to the eigensector labeled 
by (ni,ri2,P) are given by 

1*)= 2J f(xi,ai;...;x n ,a n )\xi,ai;...;x n ,a n ), 

(2) 
where the kets \xi, ot\\ . . . ; s n , a n ) denote the configura- 
tions with particles of type a» (a* = 1,2) located at the 
position Xi (Xi — 1, . . . , L) and n = n\ + n^ is the total 
number of particles. The summation {a} = {ot\, . . . , a n } 
extends over all the permutations of n numbers {1,2} in 
which n\ terms have value 1 and «2 terms the value 2, 
while the summation {x} = {x\, . . . ,x n } extends, for 
each permutation {a}, into the set of the non-decreasing 
integers satisfying x i+ \ > Xi + 1. 



The MPA asserts that the amplitudes of an arbitrary 
eigenfunction pi) can be mapped to a matrix product 



f(xi,ai; . . . ;x n ,a n ) «=^> 

J^xi - 1 ^(ai ) j^x 2 -xx- 1 j^(a 2 ) 



(3) 
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where for this map we can choose any operation on 
the matrix products that gives a non-zero c-number 
[7J. The matrices A^ are associated to the parti- 
cles of type a (a = 1,2) and the matrix E is as- 
sociated to the vacant sites [7J. Since the eigenfunc- 
tion (pk has a well defined momentum P = ^ due 
to the periodic boundary condition, the matrix product 
Pj) should satisfy the constraints f(xi, a ± ; . . . ; x n , a n ) = 
f(x\ + 1, ai;...;x n + l,a„) for x n < L, and 
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J yXi , Ot\ , . . . , X n , O-n) 
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for x„ = i. Let us consider initially the simple cases 
where n = 1 and n = 2. 

n = 1. Since the Hamiltonian (Til) is diagonal for m = 
0, let us consider only the case n\ = 1. By inserting the 
ansatz (I3J) in fl2j), the eigenvalue equation gives us 

^l.O^x-l^CD^i-x = pi 0^-2^(1)^-^1 (4) 

+r? j£; x A (1) E i - a! - 1 - e x ~ 1 a ( - 1) e l - x , 

A convenient solution of EJ is obtained by introducing 
the spectral parameter dependent matrix A^ = EA k ' 
with EA k l) = e ik A { k ] E and k e C. Inserting these rela- 
tions into Q we obtain 

e(k) = e 1 - = Tj \e-' lk + r° \e ik - 1. (5) 

The up to now free spectral parameter k is fixed by im- 
posing the boundary condition. This will be done only 
for the general n. 

n = 2. For two particles of types a,\ and ct2 {ot\, a-i — 
1,2) we have two kinds of relations coming from the 
eigenvalue equation. The configurations where the parti- 
cles are at positions (x, y) with y—x=d> 1 give us the 
generalization of Eh, and the configurations where the 
particles are at the colliding positions (y — x + 1) give us 

e ai ' a2 EA {ai) A {oi2) E = T^A^EA^E 
+T° a *lEA {ai) EA {a2) + V% *iEA (a2) A (ai) E 



the energy in terms of the spectral parameters 
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where e(k) is given by ([5]). Let us consider now ([6]) in 
the case where the particles are of same type. Using ff\ 
and pi in ffib for a\ — ct2 — 1 we obtain the relations 
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For two impurities at "colliding" positions, the eigenvalue 
equation does not fix a commutation relation between the 

(2) 

matrices A ,L . On the other hand, the sum 
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(11) 



where the number of vacant sites between the impurities 
d = y — x is a conserved charge of the Hamiltonian (nl) , 
needs to be different from zero or the MPA will produces 
an eigenfunction with null norm. Moreover, the algebraic 
expression in ( |7|) assures that any matrix product defin- 
ing our ansatz (|3]) can be expressed in terms of two single 
matrix products A^ 2) A^ 2) E L and A^A^E 1 . Using 
ufj we have, from the periodic boundary condition, 



a ( 2 ) 4( 2 ) F l 
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To satisfy this equation we should have k 2 

2irj/L (j — 0,1, ...,£— 1). Consequently, from (12) the 



most general commutation relation among the matrices 

|(2) a „j 4( 2 ) m „k v^„ m JtA 4(2) 4(2) _ 4 (2) j(2) 

1.(2) 
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_ /po qi _j_ pa 2 _|_ pai a 2 \ EA^^A^^E (6) is not zer0 ' we obtain C?' I' v 



(or S^Hkj , fc; ) = 1 in te^) by an appropriate change 

of variable in the spectral parameter k[ . Inserting these 
results in the sum of ( 11 ) and by imposing that the sum 

n 2 ) 



where we introduced T{ { = Y{{ = T^ = F^g = and 
e 0,2 = 0. A solution of these relations is obtained by 
identifying A~ a > as composed by n a spectral parameter 
dependent matrices A^' a) and A („, [13], i. e., 

A^=YeA ( H ) , with EA [a ^ =e lk > a) A {c ?lE : (7) 

with a = 1, 2 and j = 1, ...,n a . By inserting ([7]) in the 
first kind of relations, where y — x = d > 1, we obtain 



fcf^fcf 



7r(2r7i 
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= 1, 
1) 



(m = 0,l,...), (13) 



where {d v } is the set of all numbers of vacant sites be- 
tween the impurities. 

Let us consider now the case where the particles are of 
distinct species. By using ([5]), (|8| and Q, the equation 
([6]) gives us two independent relations for ai 7^ «2 = 1, 2: 
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tions need to be identically satisfied since at this level we 
want to keep k^ as free complex parameters. Then ( 14 ) 
imply special choices of the coupling constants r™ z ™: 



pi 2 _ r i o 
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with m = 0, 1, ...,L — ni — 1. By inserting ( 19 1 and using 



(10 1 in the first equation in (171 we obtain 






"LrlO , rO l p i{kj+ki) _ ik 3 

f-)" 1 - 1 ^ TT Q1 10 — 



, (20) 



The integrability conditions ( |15[ ) reduces the 2 parame- 
ters defining the Hamiltonian (111) into only 1 free param- 

pl 

eter ^§4. Moreover in order to the model be integrable, 

1 ° . 
the particle moves over impurities as if the impurities 

are vacant sites. This condition will play a fundamen- 
tal role in the spectral properties of the Hamiltonian (fil , 
giving us the full spectrum of the standard ASEP plus 



new levels. From (14) and (15) we obtain the structural 
constants: 



$i(* (2) ,* (1) ) = 
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General n. We now consider the case of arbitrary 
numbers n\, n 2 of particles of type 1 and 2. The eigen- 
value equation gives us generalizations of Q and (Rjl) . To 
solve these equations we identify the matrices Ay*) as 
composed by n a spectral dependent matrices «f\. The 
configurations where Xi+i > Xi + 1 give us the energy 
$h. The amplitudes in (pl) where a pair of particles of 
types ai and a 2 are located at the closest positions give 
us the algebraic relations Ety where the algebraic struc- 
ture constants are the diagonal 5-matrix defined by ( 10 1, 

6), S%%(kj ' ,k\ ') = 1, with coupling constants (15). 

In order to complete our solutions we should fix 



the spectral parameters (k\ 



(i) 



Using the algebraic relations Q9J) for an arbitrary ampli- 
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(17) 



where we have ( j — 1 , . . ., n%) and (q = 1 , .. ., n<z) . 

Finally, the momentum of the eigenstate is 
given by inserting the ansatz (17]) into the relation 



j \X\) ol\] . . . ; x ni OL n ) 



-ip 



f(xi + l,ai;...;x n + l,a n ): 



p = E*i 1} +E*. 
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(/ = 0,1,...,L-1). (18) 



Spectral Gap. We can write the Bethe equation ( 17) 



in a more convenient way. From the second expression 



on (17) we have 






(19) 



where j — 1, ...,ni and kj 



k 



(i) 



The Bethe equation 



( 20 1 differs from the one related to the stochastic version 
of the asymmetric XXZ chain [8] by the phase factor <f> m 
defined in (19). However, different from the phase factor 



in the asymmetric XXZ with long range interactions |10) , 
the phase <p m defined in ( 19 ) plays a fundamental role in 



the spectral properties of the model by changing the scal- 
ing behavior of the system. While in [10] the phase is a 
constant, in our model it is a function on the momen- 
tum of the state through the relation (18). For m = 
we obtain the full spectrum of the standard ASEP [S]. 
Additional energy levels are given by different choices for 
m, and the first excited state belongs to these new levels. 
Finally, it is important to notice that the Bethe equa- 



tion (|20[) do not depends on the number of impurities 
"2 



7^ 0. Consequently, the roots of (20) and the eigen- 
values ^ are independent of n 2 ^ (but the wave func- 
tion depends on n%). This huge spectrum degeneracy 
follows from the conservation of the numbers of vacant 
sites between the impurities by the Hamiltonian (nj) . Let 
us explain with the following example. Suppose we start 
with the given configuration 0120220 with one particle 
(1), 3 impurities (2) and 3 vacant sites (0). We can 
make a surjective map between all possible configura- 
tion of these particles to all possible configurations of 
a new chain with just impurities and vacant sites. For 
example 0120220 => 020220. On this new chain, we 
are looking only for the effective movement of impuri- 
ties on the chain. For simplicity, let us consider the to- 
tally asymmetric model (TASEP) where r^ \ — 1 and 
If o = 0- When the particle jumps over the impurities, 
nothing changes in the effective chain since we also have 
0210220 => 020220, then 0201220 =*- 020220, then 
0202120 => 020220, then 0202210 => 020220, then 
0202201 =*> 020220, and finally a change in the mapped 
configuration 1202200 => 202200. On other words, the 
impurities move on the mapped chain as they are just 
one "object" due to the conservation of vacant sites be- 
tween impurities. Moreover, this "object" only moves 
when the particle complete a turn over the chain. As a 
consequence, the time for the particle to complete one 
revolution is the time scale for the movements of the this 
"object" . For an arbitrary number of particles in a chain 
of length L, the time for the particles to complete one 
revolution is of order Ls. As the "object" formed by 
the impurities need to moves of order L times to span 
all possible configurations, we will need a time of order 
L'2 x L — Lz to reach the stationary state. 



The Bethe equation ( 20 1 for the TASEP can be written 



as [SJ 



(i-z 3 r(i + Zj ) 



L—n\ 



= -2 



711 1 

Li TT z l - 1 






where Zj = 2e~ lkj — 1, and we introduced an auxiliary 
complex variable Y. It can be noted that the right-hand 
side of these equations is independent of the index j, 



consequently the Bcthc equation (21 1 can be reduced to 
a simple polynomial equation of degree L as in [SJ. In 
the half-filling sector n\ — L/2 this polynomial is fur- 
ther simplified to Y = (l — z 2 ) 1 . To solve this Bethe 
equation we follow the method developed in [8J and we 
pick a Y and we take its n\th root. These roots give 
us 2ni solutions Zj. For each particular choice of n\ out 



of 2ni solutions we have to determine Y through (21). 



On the other hand, we also solve (21 ) numerically by di- 
rect computation of the roots of (21 1 up to L = 1024. 
The stationary state is provided by choosing m = and 
P = and selecting the n\ fugacities Zj with the largest 
real parts. The first excited state is also obtained by se- 
lecting the n\ fugacities with the largest real parts and 



by choosing m — 1 and P = 



1-K 



The energy gap has 



a leading behavior of KPZ L~? and a sub-leading term 
L~i related to the super-diffusion of particles 1 and sub- 
diffusion of particles 2, respectively. For m — (and also 
for n 2 = 0) the sub-leading term vanishes and we recover 
the spectrum of the ASEP without impurities 8J. For 
rn = 1 the spectral gap scales with L~i instead of L~i 
due to the sub-leading term. It is important to notice 
that the dynamics of particles 1 is not affected by the 
impurities. This explain why the model displays the full 
spectrum of the ASEP. On the other hand, the dynamics 
of the impurities is totally dependent on the particles, 
since the impurities only move when particles change po- 
sition with them. Consequently, the time to vanish the 
fluctuations on the densities of particle acts as a time 
scale for the diffusion of the impurities, resulting in a re- 
laxation time longer than the one for the standard ASEP 
and reflected in the £~i x L _1 = L~s gap. 

Eigenstates. The MPA pi) enable us to write all 
eigenstates of (nj) in a matrix product form. For a given 
solution kj and fc:- , the matrices E and -A "„, have the 

k i 

following finite-dimensional representation: 
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e ik < 
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e lk >* 
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v e lk i , 

£=1 x ' 1=1 



S\\{k h h) 
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1 



1 



1=3+1 

n 1 +n 2 

3=3+1 



where n — n\ + n 2 , I<i is the 2x2 identity matrix, and 
the dimension of the representation is 2™. The matrix 
product form of the eigenstates are given by inserting 
the matrices J7) defining the MPA ([3]) with the spectral 
parameters (20) into equation ^. As showed in [TT] . 
our MPA generalizes the steady-state Matrix Product in- 
troduced by Derrida et al [T^]. The stationary state is 



(2) 



obtained by choosing kj = and kS = in (22). How 



ever, a relation between our matrix product form for the 
steady-state and the standard Matrix Product form [T^] 
is not trivial [llj . 

Finally, like in the ASEP, the stationary state is ob- 
tained by a combination of all possible configuration of 
particles and impurities with the same probability (sat- 
isfying the order {d\, ..., d n2 } of vacant sites between im- 
purities, up to cyclic permutations). This can be seen 
from the fact that the particles move on the chain as if 
the impurities are vacant sites and the impurities moves 
on the surjective mapped chain as a single "object". 

The spectral gap for other sectors and other boundary 
conditions and the physical consequences of these, are 
currently under investigation. 

We are grateful to GM. Schiitz and F. C. Alcaraz by 
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